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We perform a systematic analysis of quantum interference in nonsequential double ionization fo¬ 
cusing on the recollision-excitation with subsequent ionization (RESI) mechanism, employing the 
strong-field approximation (SFA). We find that interference has a major influence on the shape, 
localization and symmetry of the correlated electron momentum distributions. In particular, the 
fourfold symmetry with regard to the parallel momentum components observed in previous SFA 
studies is broken. Two types of interference are observed and thoroughly analyzed, namely that 
caused by electron indistinguishability and intra-cycle events, and that stemming from different exci¬ 
tation channels. We find that interference is most prominent around the diagonal and anti-diagonal 
in the parallel-momentum plane and provide fully analytical expressions for most interference pat¬ 
terns encountered. We also show that this interference can be controlled by an appropriate choice 
of phase and excited-state geometry. This leads a to myriad of shapes for the RESI distributions 
including correlated, anti-correlated and ring-shaped. 


I. INTRODUCTION 

Quantum interference in strong-field phenomena has 
established itself as a very powerful attosecond-imaging 
tool. Well-known examples of interference are dynamic 
and structural patterns in molecular high-order harmonic 
generation (HHG) (for reviews see, e.g., [l| and our recent 
publication j2j), and the fan-shaped structure in above¬ 
threshold ionization (ATI) of rare-gas atoms 0-01- The 
connection with imaging stems from the physical mecha¬ 
nism behind these phenomena, namely the laser-induced 
recollision or recombination of an electron with its par¬ 
ent ion Q. Recombination and elastic recollision lead 
to HHG and high-order ATI, respectively. If the elec¬ 
tron is released in the continuum and reaches the detec¬ 
tor without further interaction, this will lead to direct 
ATI. Thus, HHG and ATI transition amplitudes may be 
associated with electron orbits and there will be many 
possible routes for the active electron. Hence, quantum 
mechanically, they will interfere. 

Strong-held nonsequential double ionization (NSDI) 
may also be described as the recollision of an electron 
with its parent ion. However, upon return the electron 
gives part of its kinetic energy to the core. This leads 
to the release of a second electron. Hence, a legitimate 
question is whether quantum interference also influences 
NSDI. There are, however remarkably fewer studies of in¬ 
terference effects in this context, even though NSDI has 
been investigated for over two decades. 

This apparent lack of interest may be related to the 
widespread belief that interference does not play a sig¬ 
nificant role in NSDI. This may be attributed to the suc¬ 
cess of classical NSDI models (for reviews see 0,0), 
which have reproduced key features in NDSI electron- 
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momentum distributions and have shown a very good 
agreement with experimental findings. These include, 
for instance, the shapes and maxima of the electron- 
momentum distributions, and even finer details such as 
the v-shaped structure that is a fingerprint of the long- 
range electron-electron interaction (llLll2j| . Furthermore, 
several studies of quantum-classical correspondence in 
NSDI have suggested that interference will not survive in¬ 
tegration over momentum components perpendicular to 
the laser-field polarization. This is the typical scenario in 
NSDI experiments. Such conclusions have been inferred 
from the excellent agreement between classical models 
and the full solution of the time-dependent Schrodinger 
equaHo n 0, or the strong-field approximation (SFA) 

The above-stated studies, however, have focused on 
the electron impact ionization mechanism, in which the 
first electron, by re-colliding with the core, immediately 
releases a second electron. Another physical mechanism, 
which is less studied, is re-collision-excitation with subse¬ 
quent ionization (RESI). In RESI, the first electron does 
not provide the second electron with enough energy to 
be released in the continuum. Hence, it is excited to an¬ 
other bound state, from which it subsequently leaves via 
tunneling ionization. 

The prevalent view is that the above-mentioned time 
delay leads to anti-correlated electron momentum distri¬ 
butions. They would populate the second and fourth 
quadrant of the plane spanned by the electron momen¬ 
tum components p n \\, (n = 1,2) parallel to the driving- 
field polarization. These features have been identified in 
experiments performed in the below-threshold intensity 
regime |17H2ll |. However, more recent results strongly 
suggest that this interpretation is over-simplified. 

For instance, recent experiments of RESI of Ar with 
few-cycle pulses have revealed cross-shaped distributions 
strongly localized along the axes p n \\ = 0 0- Subse¬ 
quent studies by the same group have shown that if the 
pulse length is increased, the distributions spread across 
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the four quadrants of the parallel momentum plane with 
a slight preference for back-to-back emission [ 23 . This 
agrees with the findings in [HI for NSDI in Xe, namely 
distributions equally occupying all momentum quadrants 
that exhibited RESI characteristics. Hence, it is plausible 
that in fl8l - 0 , the signal in the first and third quadrant 
of the parallel momentum plane, which was dismissed as 
electron-impact ionization, could in fact be RESI. 

This affirmative is backed by theoretical findings using 
methods as diverse as the SFA and related methods |25l - 
0 and classical-trajectory f2il. [28i [dp| computations. Al¬ 
though back-to-back emission was highlighted, in many 
classical-trajectory studies cross- or ring-shaped distribu¬ 
tio ns sp reading across all quadrants have been identified 
[24 !28l - [30[ . This behavior has even been found for inten¬ 
sities far below the threshold, for which electron-impact 
ionization can definitely be ruled out [0. In particular 
for semi-analytical methods such as the SFA 0, [26j, 0 
or the quantitative rescattering theory 0, exclusively 
fourfold symmetry distributions were identified for RESI. 
One should note, however, that in none of the SFA com¬ 
putations has interference between different events or ex¬ 
citation channels been incorporated. Recently, however, 
SFA computations using inter-channel interference have 
shown that the fourfold symmetry can be broken ET2. 
Indeed, it has been argued that quantum interference is 
paramount, and that it may lead to anti-correlated dis¬ 
tributions. 

Nonetheless, more systematic studies of quantum- 
interference effects in RESI are missing. This interfer¬ 
ence may occur (a) between events which are displaced 
by half a cycle and those present due to the symmetry of 
indistinguishable electrons, (b) between different chan¬ 
nels of excitation for the second electron and (c) between 
different orbits along which the first electron may return. 

In this article, we present a systematic analysis of 
the two first types of interference, employing the SFA. 
The SFA, if used in conjunction with the steepest de¬ 
scent method, provides a very intuitive interpretation in 
terms of electron orbits, and retains quantum interfer¬ 
ence and tunneling. This makes it an ideal tool for an¬ 
alyzing different types of interference. The last type of 
interference is incorporated in the model, but is washed 
out when the perpendicular momentum components are 
integrated over. This has been studied in previous pub¬ 
lications 0. IT >j. 

This work is organized as follows. In Sec. ED we briefly 
review the necessary background for understanding the 
subsequent results. These are provided in Secs. ED El 
and 0 In Sec. Em we investigate symmetry-related in¬ 
terference and provide expressions for the features en¬ 
countered. Subsequently, in Sec. IIV1 we analyze how the 
geometry of the bound states involved modifies this in¬ 
terference, and in Sec. 0 we study how different excita¬ 
tion channels interfere and how this affects the electron- 
momentum distributions. Finally, in Sec. IVI1 we provide 
an overall discussion and state our main conclusions. 


II. BACKGROUND 

A. Transition amplitude and saddle-point 
equations 

Within the strong-field approximation and in atomic 
units, the RESI transition amplitude reads 


M(pi,p 2 ) = dt dt dt d 3 k 
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incorporate all information about the interactions and 
electronic bound states. Specifically, Eqs. ©, © and 
© are related to, the ionization of the first electron, the 
recollision of the first electron with excitation of the sec¬ 
ond electron, and the tunnel ionization of the second elec¬ 
tron, respectively. Therein, V(ri) and Vion^) denotes 
the binding potential “seen” by the first and the second 
electron, respectively, and V'i 2 (i'i,r 2 ) gives the electron- 
electron interaction. Furthermore, k(r) = k + A(r) and 
p n (r) = p n + A(r) (t = t, t',t") in the length gauge, and 
k(r) = k and p„(r) = p„ in the velocity gauge, with 
n = 1,2. In our previous publication (2f|, we have veri¬ 
fied that, in practice, the results obtained in both gauges 
lead to qualitatively similar results. Here we employ the 
latter gauge in order to avoid bound-state singularities. 

The transition amplitude © describes a process in 
which the first electron, initially bound in the ground 
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state \ip[ 9 ^), is released at a time t" into a continuum 
state, which is approximated by the Volkov state |k(f")). 
Subsequently, it remains in the continuum from the time 
t" to the time t' with intermediate momentum k. At t', 
it returns to its parent ion and interacts with a core elec¬ 
tron via Vi 2 * This interaction excites the second electron 
from the ground state | ip 2 3 ^) °f the singly ionized target 
to the state \ip 2 e ^)- The first electron reaches the detector 
with final momentum pi immediately after rescattering. 
The second electron remains bound until a later time t, 
when it is released by tunnel ionization into a Volkov 
state |p 2 (f)). It reaches the detector with final momen¬ 
tum p 2 . The ground-state energy of the neutral system 
is given by Ei g , and the energies of the ground and ex¬ 
cited states of the singly ionized target are E 2g and E 2e , 
respectively. For details on how this transition amplitude 
is derived we refer to our previous publication |2fj . 

Throughout, we employ the steepest descent method 
in order to compute the transition amplitude ©■ In this 
method, we seek variables t, t't" and k so that the action 
is stationary. This leads to the saddle-point equations 


[k + A(f")] 2 = —2£i g , 

(6) 

i r*' 

k t'-t" j t „ dTA{T) ’ 

(7) 

[pi + A(t') 2 ] = [k + A(t')] 2 - 2 (E 2g - E 2e ), 

(8) 

and 

[p 2 +A(f)] 2 = A2£ 2 e. 

(9) 


Eqs. © and © give the energy conservation of the first 
and second electron at the instant of tunnel ionization. 
For the former, this occurs from the ground state at the 
time t", while the latter tunnels from an excited state 
at a later time t. Neither has a real solution, which re¬ 
flects the fact that tunneling has no classical counterpart. 
Eq. 0 restricts the intermediate momentum of the first 
electron so that it returns to the core, which is assumed 
to be located at the origin. Finally, Eq. © states that, 
upon return, the first electron gives part of its kinetic 
energy upon return to “bridge” the gap E 2g — E 2e and 
promotes the second electron to an excited state. We 
use both the standard saddle point approximation and 
a uniform asymptotic expansion whose only applicability 
requirement is that the orbits occur in pairs. For details 
on these methods see our previous work [ 34 } . 

B. Momentum constraints 

For simplicity, we will consider a monochromatic field 
throughout. Hence, the vector potential reads 

A(t) = 2 ^/Up cos (ujT)e\\, (10) 


where r is a generic time that may be t,t' or t" , u is 
the driving-field frequency and U p = //(4w 2 ) is the pon- 
deromotive energy, which is proportional to the intensity 
I of the driving field. This choice of field means that 
A(t ± T/ 2) = —A(t), where T = 2'k/uj denotes the field 
cycle. Thus, events whose times are displaced by half a 
cycle are related by momentum inversion. Furthermore, 
since both electrons are indistinguishable, one must also 
exchange pi and P2 in Eq. © and add the corresponding 
amplitudes. 

The distributions will be located around (pi,p 2 ) = 
(±2y/f7p, 0). This comes from the fact that the rescatter¬ 
ing of the first electron and ionization of the second elec¬ 
tron occur most probably near field crossings and crests, 
respectively. 

Estimates for the regions in the parallel momentum 
plane to be populated follow from the saddle-point equa¬ 
tions. For the first electron, Eq. © gives 

±2y/%- v / 2A^<pi|, <±2 ^+v / 2A^, (11) 

where A E = E]^ n {t', t") — E exc yields the energy differ¬ 
ence between the kinetic energy E kin (t', t") of the first 
electron upon return and the energy E exc = E 2g — E 2e + 
pl ± / 2 . The above-stated inequality indicates that the re¬ 
gion where rescattering has a classical counterpart, which 
is largest if pn = 0. For the second electron, one must 
bear in mind that Eq. © is formally identical to that de¬ 
scribing tunnel ionization in direct ATI. The direct ATI 
cutoff energy is 2 U p , so that 

— 2 \/Up ^ P2|| ^ 2 \/V p . (12) 

In this latter estimate, we have considered p 2 j_ = 0. 
Eqs. (fill) and (fl2l) give cross-shaped electron momen¬ 
tum distributions strongly located around the axes of the 
Pi\\p 2 \\ plane. Detailed explanations of these constraints 
have been provided elsewhere [25], H(| [35} . 

III. INTERFERENCE OF EVENTS 

Here we will analyze interference between events dis¬ 
placed by half a cycle and those present due to the 
particle exchange symmetry of the system. This leads 
to four transition amplitudes, M(pi,p 2 ), M(p 2 ,pi), 
M(— Pi,—p 2 ) and M(— p 2 ,—Pi), which must be com¬ 
bined. Due to the localization of these transition ampli¬ 
tudes near the negative half axis, positive p 2 || half 
axis, positive piy and negative p 2 || half axis, i.e., occupy¬ 
ing the left, upper, right and lower regions in the parallel 
momentum plane, we relabel them M;, M u , M r and 
respectively. A schematic representation of the momen¬ 
tum regions occupied by the different transition ampli¬ 
tudes is provided in Fig. |T] 

Throughout this analysis we will compare coherent 
and incoherent sums of these amplitudes integrated over 
momentum components perpendicular to the laser held, 
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FIG. 1: Schematic representation of the momentum regions 
occupied by the transition amplitudes Mi, M u , Md and M r , 
which are displayed as the rectangular patterned regions. The 
overlap regions indicate areas within these constraints for 
which quantum interference may occur. The intensity repre¬ 
sented in the figure is high enough to allow some interference 
at the origin (pi||,p 2 ||) = (0,0). 


which are given by 

W (pi|| ,P2\\) = Jd 2 pi±d 2 p2±\Mi + M u + M r + Md \ 2 (13) 

and 

W(pn\,p 2 \i) =/ d 2 p 1± d 2 p 2 ± 

x (|M,| 2 + \M U \ 2 + \M r \ 2 + \M d \ 2 ), (14) 

respectively. 

Quantum interference occurs predominantly in the 
overlap regions in Fig. [I] Apart from the region around 
(pi ||5 P 2 \\ ) = (0,0), in which, potentially, all amplitudes 
may interfere, due to the constraints discussed in the pre¬ 
vious section, we expect that Mi and M u will interfere 
predominantly in the second quadrant and that Mi and 
M d will overlap in the third quadrant of the Pi\\P 2 \\ plane. 
Similarly, interference between M r and M u is expected to 
occur in the first quadrant, and between M r and Mi will 
take place mostly in the fourth quadrant. For simplicity, 
throughout this section we will neglect the prefactors in 
Eqs. (1131) and m- This will help us identify how the 
phases determined by the corresponding actions interact 
without further momentum bias. The prefactors will be 
reintroduced in Sec. m 

In Fig. [2j we display the coherent and incoherent sum 
for three driving-field intensities. The figure shows that 
the interference between different events survives the in¬ 
tegration over the transverse momentum coordinates, as 



FIG. 2: Coherent and incoherent sums of all amplitudes inte¬ 
grated over the perpendicular momenta as given by Eq. USD 
and CD- The columns from left to right show a coherent sum 
of amplitudes [panels (a), (d) and (g)], an incoherent sum of 
amplitudes [panels (b), (e) and (h)], and the difference be¬ 
tween the two [panels (c), (f) and (i)[. The rows show differ¬ 
ent laser intensities of I = 2.28, 4.56 and 6.84 (x 10 13 W/cm 2 ) 
from top to bottom with values for U p of 0.05, 0.1 and 0.15 
a.u. corresponding to an angular frequency u> = 0.057 a.u. 
The RESI channel corresponds to a ground state valence or¬ 
bital 3s and an excited orbital 3p for the second electrons. 
The ionization potentials of the ground state of the first and 
second electron and the excited state of the second electron 
are I\ g = 0.58, h g = 1.02 and / 2e = 0.52 a.u. respectively. 
The diagonal and anti-diagonals piy = ±p 2 [| are indicated 
with the orange lines in the figure. The signal in each panel 
has been normalized with regard to its maximum. 


there are obvious differences between coherent and in¬ 
coherent sums of events. Clearer features can be out¬ 
lined from the difference of the two probability maps. 
There are maxima along the diagonal and anti-diagonal 
at all intensities and hyperbolic fringes whose presence 
becomes more obvious as the intensity increases. For 
higher intensities the patterns become more complicated. 

These features can be explained by looking at the in¬ 
tegrand of the coherent sum, which can be rewritten in 
terms of the actions Si, S r , S u and S d associated with 
the above-stated amplitudes. A common factor can be 
taken out, leaving terms that will contribute to the in¬ 
terference. Explicitly, 


ft(Pi\\,P2\\)= / d A p± 


I d 3 t(e 


iSl _J_ g*Sr 


JSu 


„iSd'\ 


(15) 



<?te iS '( 1 + e lo ‘ lr + e*“ ! “ + e ia ‘ d ) 


where the action Si = S(pi, p 2 , k, t, t 1 , t") is associated 
with the matrix element Mi = M(pi, P 2 ) giving the left 
peak. The integrals over time and momenta have been 
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abbreviated as 

/ /•OO nt nt' 

d 3 t = dt dt dt , (16) 

J— oo J — oo J—oo 

and 

J d A p±= J J d 2 p 1 ±d 2 p 2 ±, (17) 

and the phase differences between the actions read 
a-u = \{p\~ pl) (t - t') 

+ 2 ^ ,i|| -P 2 ||)(sin(urt) - sin(wt')), (18) 

aiu = + 2E2e + 2Ei g + Pi + P 2 ) 

+ \ (p? - pI) (t ~ *0 

- 2 ^ P (Pi\\ +P 2 ||)(sin(wt) - sin(wt')) (19) 

UJ 

and 

air = — (4 Up + 2E2 e + 2Ei s + Pi + P 2 ) 

ZUJ 

-^ (pi|| sin(wf') +P 2 \\ sin(wt)) . ( 20 ) 

UJ 

In Eq. (fl8l) - (l20l) . the general form lias been 
adopted, where the indices i,j refer to the interfering 
amplitudes. Constructive interference requires that the 
integrand in Eq. m is maximized, which will occur 
when onj = 0, or as small as possible. We will start by 
investigating the “left-down” phase difference (fl 8 l) . be¬ 
tween the actions associated with the left and lower peak. 
This phase vanishes for arbitrary times t, t! if conditions 

P?-P 2=0 (21) and Pi\\ = p 2 \\ (22) 


are satisfied. 

Condition m, if written as a function of the paral¬ 
lel and perpendicular momentum components for p 2 ± ~ 
Pi_l ^ 0 , give the hyperbolae 


Pill _ P2|| 

pIj_-Pi± P 2 2± - Pl± 


(23) 


whose asymptotes lie at the diagonal and anti-diagonal 
Pi || = ±P 2 || and whose vertices and transverse axis will 
depend on whether p^ ± — pf ± are positive or negative. 
The former and latter case will lead to hyperbolae with 
transverse axes along pi\\ andp 2 ||i respectively. For equal 
transverse momenta, instead, condition m will give 
Pi \| = ±p 2 || 1 he., the diagonal and the anti-diagonal. 
In this case, the interference condition is independent 
of the transverse momenta, so that they are expected to 



FIG. 3: Correlated electron-momentum distributions ob¬ 
tained by combining the transition amplitudes Mi and Md, 
isolating the effect of aid, integrated over the transverse- 
momentum components. The left, middle and right columns 
have been computed for laser intensities of I = 2.28, 4.56 and 
6.84 (xlO 13 W/cm 2 ), with values for U p of 0.05, 0.1 and 0.15 
a.u., respectively. Panels (a), (b) and (c) show the coherent 
sum flcoh = \Mi + Md | 2 , while in panels (d), (e) and (f) the 
incoherent sum = \Mi\ 2 + |JWd| 2 is displayed. The in¬ 
tensities and ionization potentials are the same as in Fig. [2] 
The diagonals 'pu\ = p 2 || are indicated with the orange lines 
in the figure. The signal in each panel has been normalized 
with regard to its maximum. 


survive when the integration over these variables is per¬ 
formed. The hyperbolae, on the other hand, depend on 
the transverse momentum coordinates, but may survive 
integration. If this happens, however, integration may 
influence their transverse axes, vertices and foci. 

The analysis performed above suggests that there will 
be maxima along the diagonal and the anti-diagonal, and 
that there could be hyperbolic fringes in the coherent sum 
of the two-electron transition amplitudes, in agreement 
with Fig. [2j In Fig. [3l we have a closer look at this inter¬ 
ference, and plot a partial distribution in which only Mi 
and Md summed, coherently and incoherently (upper and 
lower panels, respectively). The strongest feature in the 
figure is the maximum along the diagonal, which comes 
from condition (l22l) and also from the case pu = p 2 ±. 
related to the hyperbolic condition (1211) . Parallel to the 
diagonal, there are also interference fringes, whose num¬ 
ber increases with the driving-field intensity. The inter¬ 
ference maxima along the anti-diagonal cannot be seen 
as the partial sum employed in the figure is vanishingly 
small in the second quadrant of the parallel momentum 
plane. 

An estimate for the position of the fringes can be ob¬ 
tained by considering the coherent superposition of Mi 
and Md, and expanding the momenta in the vicinity of 
the diagonal, i.e., pi y = — p 2 \\ + S. Fringes will occur for 
exp [iaid] = ± 1 , i.e., for aid = mr, where even and odd n 
give maxima and minima, respectively. Assuming small 
momenta p„(n = 1 , 2 ), rescattering times at field cross¬ 
ings and ionization times at the subsequent field crest 
[■ t' = mr/uj and t = (2 n+ 1)7t/(2w)], the fringe spacing 
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FIG. 4: Absolute value of the difference between the sums 
from the upper and lower panels in Fig. [3] Panels (a) and (b) 
correspond to U v = 0.05 and U p = 0.1 a.u respectively. The 
approximate fringes given by Eq. (12411 are marked by black 
lines. The panels have been plotted in a logarithmic scale 

(iog(|n£ h - ntfi)). 


can be approximated as 


\S\ 


UJT17T 


(24) 


The above-stated equation shows that the spacing be¬ 
tween the fringes is inversely proportional to the driving- 
field strength, in agreement with what has been observed 
in the previous figures. 

The interference patterns are highlighted in Fig. |4j 
where we display the difference between the coherent and 
the incoherent sum, for the two lower driving-field inten¬ 
sities in the previous figure. Overall, for small momenta 
the fringe spacing exhibits a very good agreement with 
Eq. (l24l) . Furthermore, all panels in the figure exhibit 
clear hyperbolic structures, whose number increases with 
the driving-field intensity. One should note that their 
transverse axis is not located along p n n = 0 (n = 1,2). 
This displacement is probably related to the integration 
over the transverse momenta, which influence the direc¬ 
tion of the hyperbolae. Furthermore, the last diagonal 
term will act to shift the center of the hyperbola along 
the diagonal, which can be observed by the fact the hy¬ 
perbolae are opening, instead of exhibiting asymptotic 
behavior towards the diagonals. As the laser intensity 
increases the hyperbola should be shifted further from 
(pi||,p 2 ||) = (0, 0) and the number of fringes increases, as 
indicated by Eq. (liMl) . 

The remaining phase shifts, cq r and aid, will not van¬ 
ish. However, by an adequate choice of parameters one 
may identify momentum regions in which they are small¬ 
est, which will give rise to interference maxima. The ex¬ 
ponent oliu , which gives the interference between the left 
and the upper peaks, behaves in a similar way as aid, 
with the main difference that Eq. m shows an addi¬ 
tional phase, with regard to Eq. (TtSl) giving the left-down 
phase difference. This phase depends on pi +p| and has 
a constant factor. Furthermore, the last term in Eq. m 
causes a strong enhancement along the anti-diagonal. 


2 

CN 

Q.1 

qT-i 

-2 

2 

^ 0.1 

qT-1 
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FIG. 5: Correlated electron-momentum distributions ob¬ 
tained by combining the transition amplitudes Mi and M u , 
isolating the effect of ai u , integrated over the transverse- 
momentum components. The driving-field parameters and 
ionization potentials are the same as in Fig. [3] We have also 
employed the same normalization and labeling as in Fig. [3] 
with the coherent and incoherent sums in the upper and lower 
panels, respectively. The anti-diagonals pin = — p 2 || are indi¬ 
cated with the orange lines in the figure. 



The coherent and incoherent sums of Mi and M u are 
presented in Fig. [5] The figure shows a very clear max¬ 
imum along the anti-diagonal = — p 2 \\, and interfer¬ 
ence fringes with a richer substructure than the previous 
partial map. These effects are caused by the additional 
phases mentioned above. An estimate of the position of 
the fringes is not straightforward. However, we have ver¬ 
ified that their spacing, for small momenta, is approx¬ 
imately one fourth of that given by Eq. (l24l) . It also 
decreases with driving-field intensity. 

Finally, we display the partial sum between the left 
and right amplitudes M r and Mi (Fig. [5]). In this case, 
the interference effects are minimal and only present close 
to the origin (pi||,P 2 ||) = (0,0). This is expected as ac¬ 
cording to the constraints the overlap between both am¬ 
plitudes is vanishing in the other momentum regions. In 
the case of cq r there is no obvious condition on the paral¬ 
lel momentum, independent of time or the perpendicular 
components, other than pu i and p 2 \\ being close or equal 
to zero. A time-dependent condition can be extracted, 


sin^t') 

V2\\ =- . / ,N -Pl||- 

sm(wf) " 


(25) 


These trajectories will overlap for very low values of par¬ 
allel momenta. In this case, for the dominant trajectories 
t' is near a crossing and t is near the next maximum [2(1] . 
These trajectories are located close to the axis. Near a 
crossing, sin(wt / ) ~ Lot' and near a maximum sin(wt) ~ 1. 
This strongly suggests that the slope in the overlap re¬ 
gion will be constant as the rescattering time will not 
vary substantially. This is approximately the behavior 
observed in Fig. [G] 
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FIG. 6: Correlated electron-momentum distributions ob¬ 
tained by combining the transition amplitudes Mi and M r , 
isolating the effect of a; r and integrating over the transverse 
momentum. The field parameters, ionization potentials and 
plotting style are the same as in Fig. [3] 


Channel Excited-State Configuration E 2 e (a.u.) 


1 

3s3p 6 (3s —> 3 p ) 

0.52 

2 

3p 5 3d (3p —> 3d) 

0.41 

3 

3p*4s (3p —> 4s) 

0.4 

4 

3p s 4p (3p —> 4p) 

0.31 

5 

3p 5 4d (3p 4d) 

0.18 

6 

3p 5 5s (3p —± 5s) 

0.19 


TABLE I: Relevant excitation channels for Ar + , ordered ac¬ 
cording to principal and orbital quantum numbers for the sec¬ 
ond electron’s excited state. From left to right, the columns 
give the number associated with the channel, the electronic 
configurations for the sub-levels involved in the excitation 
and the absolute value E 2e of the excited-state energy, re¬ 
spectively. For clarity, the excitation pathway for the second 
electron is given in brackets. 


IV. THE EFFECT OF THE PREFACTORS 


Additionally to the interference effects studied above, 
the prefactors £3} and m will introduce a momentum 
bias, which influences the shapes and, in principle, the 
quantum interference between events or channels. In the 
specific problem addressed in ]32j, the target chosen is 
Argon, whose first and second ionization potentials are 
Ei g = 0.58 a.u. and E 2g = 1.02 a.u., respectively. For 
the parameter range of interest, there exist six relevant 
excitation channels, which are provided in Table [T] and 
involve excitations to states of very different spatial ge¬ 
ometry. Hence, they will give us a fairly good idea about 
the role of the prefactors. Throughout, we will restrict 
our studies to m = 0, in order to facilitate a comparison 
with the results in 13211. 


In Eq. (l26l) and (l27l) . we give the general expressions 
for the excitation and ionization prefactors, respectively, 
assuming that all bound-states involved are of the form 
ipnim( r) = i.e., hydrogenic states. These 

prefactors have been first derived in [251 ] in the context of 
a qualitative analysis, so that only their functional form 
has been emphasized. In the present work, we go beyond 
those qualitative expressions and include all normaliza¬ 
tion constants and phases, as they will be necessary for 
computing coherent superpositions. 
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Now here is the expression for V P2e again with all nor¬ 
malization constants and phases. 
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(27) 

The above-stated prefactors have radial and angular 
nodes. For the first electron, Eq. (ESI) depends on the 
intermediate momentum k(t'',t'), which will vary with 
regard to pi. This will lead to these nodes being washed 
out to a great extent. We have verified that this happens 
even if the integration over pu_ is not performed. In gen¬ 
eral, transverse momentum integration will cause further 
blurring. Mostly, the prefactor Vp ie .kg will cause a shift 
in the peaks of the electron momentum distribution from 
PiII = zL2y/Up and alter their width. 

The effects caused by the prefactor V P2e are much more 
dramatic. This has been observed in our previous pub¬ 
lications [Iji, [36j for atoms and molecules, but has not 
been investigated systematically. Similarly to what is 
observed for hydrogenic wave functions, the number of 
radial nodes is given by n e — l e — 1, and angular nodes by 
l e . This is because, formally, the prefactor is the Fourier 
transform of a hydrogenic excited state ipnj e m e (r 2 ) mod¬ 
ified by the interaction I^ 0 n(r 2 ) = 1/^2- Since Vi on and 
exp(*p 2 ) ■ r 2 have no nodes, the number of nodes will 
be preserved but their energy positions will be differ¬ 
ent, if compared to the momentum-space wave function 
'tpn e i e m e (p 2 )- Their number and position with regard to 
the momentum p 2 = \ p\\\ +p!± are given i n table [TT] 

According to tabic |TT1 the radial nodes will manifest 
themselves as circles in the P 2 ||P 2 _l plane. They are 
clearly seen if we fix the momentum of the first electron at 
(Pi\\ iPi±) = (2y/I7p, 0) and plot the probability distribu¬ 
tion as a function of the momentum components P 2 || and 
P 2 _l of the second electron. This procedure is similar to 
the computation of partial momentum maps employed in 
our previous publications [33| [35[, and provide a wealth 
of detail which is lost if the transverse momentum inte¬ 
gration is performed. 

Fig- 0 displays these distributions for the six chan¬ 
nels in Tabic Q] The panel labels each correspond to the 
channel number, which is detailed in table |U The cir¬ 
cle P 2 || + P 2 ± = 4 Up indicates the direct ATI cutoff, ac¬ 
cording to the condition (fTEI) . Changes in the shapes 
of the distributions will be caused by nodes within this 


FIG. 7: Cross-section of the total probability distribution 
with pi fixed at (pi||,Pi_l) = (2y / Cp, 0), which gives an effec¬ 
tive partial probability distribution over P 2 . The ponderomo- 
tive energy is given by U p = 0.1 a.u.(/ = 4.56 x 10 13 W/cm 2 ). 
A logarithmic scale has been used to highlight the orbital- 
geometrical features. The radial and angular nodes result¬ 
ing from the second ionisation prefactor are marked by green 
circles and red line respectively. The direct ATI cutoff 
P 2 || + P 2 ± = 417p is marked with a dashed circle. Beyond 
this point the probability distribution decays exponentially. 
Phases for each prefactor are indicated by + and - signs, with 
a change in sign indicating a flip. The signal in each panel 
has been normalized with regard to its maximum value. 


region. The radial nodes will then be particularly im¬ 
portant for highly excited states, as in this case E^e is 
small. Physically, this is related to the fact that localiza¬ 
tion in momentum space corresponds to a position-space 
spread. The smaller the binding energy, the more delo¬ 
calized i>nj e m e (^ 2 ) will be. 

The effect of the radial nodes can be seen by compar¬ 
ing channels 3 and 6, which involve s states. For channel 
3, there is only one radial node in the momentum region 
of interest, while for channel 6 the two existing nodes in¬ 
fluence the electron momentum distributions. This will 
lead to an overall narrowing in momentum space. In 
the remaining channels, additionally to this effect, there 
are also angular nodes, which behave in very distinct 
ways. For l e = 1 (channels 1 and 4) they lead to a 
strong suppression in the electron-momentum distribu¬ 
tions for p 2 || = 0. Since these nodes occur for all P 2 ±, 
they will survive the transverse-momentum integration. 
This will cause the correlated two-electron distributions 
to move away from the axes. For d states, there are x- 
shaped nodes which intersect at {jp 2 \\iPi\\) = (0,0). We 
have verified that these nodes will also survive the inte¬ 
gration over the transverse momentum components, but 
will lead to a secondary, much weaker maximum at the 
axes instead of a complete suppression. 

In Fig. [SI we plot the incoherently symmetrized, cor¬ 
related distributions, for the same channels as in Fig. 0 
The figure in fact shows an overwhelming influence of 
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Channel 

and 

State 

^71 e l e 

Numerator 

Roots 

1 3 p 

1 

2E2e ~ V2 

P2 = \J2E2r 

2 3 d 

0 

const. 

no roots 

3 4s 

3 

8 E\ £ - 28E\ B p\ + 14 E 2e p 4 2 - p\ 

P2 = V2E2e, P2 = \j (/> ± Ay/2)E 2e 

4 4 p 

2 

20E% e - 28 E 2e p% + 5 pj 

P2 = ^2 (7 ± 2\fE) E 2e 

5 4 d 

1 

2-E'2e — v\ 

P2 = \/2E2e 

6 5s 

4 

80 Et - 480Ef e p| + 504E| e pi - 120 E 2e p 6 2 + 5pl 

p 2 = J (2 ± E 2e , P2 = \J (10 ± Ay/E) E 2e 


TABLE II: Number of radial nodes b ne i B = n e — l e — 1 and the numerator polynomials (and their associated roots) that give 
rise to these nodes. Note that pi = p\ || + P 2 _l, so that the expressions in the third column describe circles in the P 211 P 2 X plane. 
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FIG. 8: Full probability distribution with both prefactors 
included, an incoherent sum of events has been used. The 
panels are marked with the channel number and excitation 
state (top left and right corners, respectively). Red lines 
have been used to show the splitting caused by angular nodes 
and green circles mark the secondary peaks due to radial 
nodes. The ponderomotive energy is given by U p = 0.1 
a.u.(J = 4.56 x 10 13 W/cm 2 ). The yield in channel four has 
been over-exposed in order to show the secondary nodes. 


the prefactor V P2e . Angular nodes in p 2 \\ are clearly vis¬ 
ible as cuts marked by orange lines and radial nodes can 
be seen by small secondary peaks marked by green cir¬ 
cles. Only for very loosely bound states does the exci¬ 
tation prefactor lead to some substructure (see channel 
4), although it is an order of magnitude below the main 
peak. This figure establishes which substructure comes 
from the prefactors themselves, so that they cannot be 
attributed to the interference between different events. 

If a coherent sum is considered upon symmetrization, 
Fig. [9l the same diagonal fringes can be seen as in 
Fig. [2] Given interference only occurs along the diago¬ 
nals, localization for s states by the two prefactors, which 
narrows the distribution width and pushes the peak away 
from the origin, cause the diagonal and central region to 
be minimally occupied. Hence, little interference occurs 


FIG. 9: Full probability distribution with both prefactors in¬ 
cluded, for identical parameters as in figure [8l except that a 
coherent sum of events has been used. The same panel la¬ 
bels are used as in Fig. [8] The signal in each panel has been 
normalized with regard to its maximum value. 

for s states. For p and d states, there is a lot of inter¬ 
ference as the effect of the angular nodes is to split the 
distribution apart, widening it causing much of it to be 
along the diagonal. The actual type of interference is 
unchanged from figure [2j we verified this by looking at 
the phase information from the prefactor. For the sec¬ 
ond ionization prefactor, looking a the phase plotted over 
p 2 , the nodes represent a phase shifts of 7r, if this is ap¬ 
plied to the partial momentum distribution there is little 
change in the resulting phase map and the effect of this 
after integration over p± will be lost entirely. Hence the 
prefactors effect the interference only by localization and 
all the effects derived, discussed earlier, are still valid. 


V. INTERFERENCE OF CHANNELS 

We will now study the quantum interference between 
the different excitation channels in Table |U A uniform 
superposition of channels is used, which can be justified 
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FIG. 11: Channel sum 1 and 4 without prefactors. Panel 
a) shows an incoherent sum of channel and a coherent sum of 
events, panel b) shows a coherent sum of channels and events, 
panel c) is the difference between the two. The driving-field 
parameters are the same in Fig. [8] The signal in each panel 
has been normalized to its maximum value. 


FIG. 10: Full coherent and incoherent superpositions of all 
channels in TableUJ for the same field parameters as in Fig. [5] 
Panel (a) is an incoherent sum of all the channels and their 
events, panel (b) is an incoherent sum of channels with a co¬ 
herent sum of event, panel (c) is a coherent sum of channels 
with an incoherent sum of events and panel (d) is a coher¬ 
ent sum of channels and events. The symbols i and c in the 
bottom right corners denote incoherent and coherent sums for 
event and channel, respectively, with event preceding channel. 
All panels use the same arbitrary scale. 


if one views each channel as a path the second electron 
can take from its ground state to the final Volkov state. 
Hence, the final transition amplitude should sum over the 
possible channels, leading to | J0 C M C | 2 , where M c is the 
transition amplitude calculated for each channel. 

In Fig.[Tn]we plot the full sum of channels 1 to 6 , using 
different combinations of coherent and incoherent super¬ 
positions for events and channels. The figure shows that a 
fourfold momentum symmetry only occurs if the channels 
and events are summed incoherently [panel (a)]. Once 
quantum interference is introduced, only the reflection 
symmetry with regard to the diagonal or anti-diagonal 
remains, as shown in panels (b) to (d). In this case, the 
features along the diagonal and the anti-diagonal differ. 
However, only channel interference [panel (c)] exhibits 
a diagonal enhancement. The anti-diagonal fringes only 
come from event interference [see panel (b)]. The diago¬ 
nal enhancement and breaking of symmetry in panel (c) 
is consistent with what was found in j32j . 

A legitimate question is whether one may identify dom¬ 
inant channels and/or features related to the channel 
type in the superpositions presented above. The shapes 
of the superpositions in [TO] suggest that excitations in¬ 
volving p and d states prevail. Table Hill shows the mean 
values of the correlated electron momentum distributions 
for each channel, which are comparable. Since one chan¬ 
nel does not dominate significantly over the rest, inter¬ 
ference is expected to be important. This is contrary 
to the results in [32|, where channels 1-3 were found to 
dominate. 

More insight is obtained by considering superpositions 


of two channels, which may be incoherent or coherent. 
The former and the latter case are given by 

^/n(Pi||,P2||) = l|M 1 || 2 + ||M 2 || 2 (28) 

and 

^Coh(Pi\\,P 2 \\) = \\M\ +e^M 2 || 2 , (29) 

respectively. In the coherent sum (l29l) . we have included 
a phase (f) that can be used to manipulate interference 
effects such as diagonal or anti-diagonal enhancement. 

Without the effect of prefactors there is little qualita¬ 
tive difference between the possible channel sums, given 
that the actions only differ by the term E 2 g t [see pan¬ 
els (a) and (b) in Fig. QTj. Nonetheless, in the difference 
between the coherent and incoherent sums we can see hy¬ 
perbolic fringes [Fig.QjJc)]. We have verified empirically 
that the position of the fringes is determined by the value 
of the phase 4> and the thickness of the fringes is inversely 
related to _E 2ei — I? 2e2 . Due to their location in the P 2 ||Pi|| 
plane, the most significant interfering terms will be equiv¬ 
alent events between channels, e.g., Mu and M 2 ;, but the 
terms related by particle exchange such as Mu and M^d 
will also be important. The prefactors add a momentum- 
dependent phase difference between the two channels. In 
the case of the second ionization prefactor, which mainly 
determines the interference effects, the phase is constant 
but inverts when a nodal line is crossed. The other pref¬ 
actors depend on k, which has a complex phase relation 
determined by the saddle point equations. 

Figure fl2l shows a selection of particular interference 
phenomena occurring in two-channel sums. In panel (a) 
the recollision and second ionization prefactor both cause 
phase inversion, which overall cancels. The diagonal en¬ 
hancement comes from the fringes related to channel in¬ 
terference shifted by a phase of —i, whereas i would cause 
an anti-diagonal enhancement. For panel (b) the total 
effect of the prefactors is to cause an inversion in the 
channel-interference fringes. This leads to a suppressed 
signal along the diagonal. However, the effect is not as 
strong as the fringes are distorted by V P2e . In both pan¬ 
els (a) and (b), the thickness of the channel-interference 
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1 
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3 
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3 

4.94 x 10 -23 

3 

2.52 x 10“ 2 ° 

6 

6.17 x 10 _3 ° 

6 

1.21 x 10 -23 

6 
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TABLE III: Mean values of the two electron parallel momentum probability distribution of each channel for different laser 
intensities, within the parameter range of interest. These distributions have been computed for a monochromatic field. 



FIG. 12: Two-channel sums with prefactors, for the same 
driving-field parameters in Fig. [8] The numbers at the top 
right in each panel labels the excited states used in the super¬ 
position, with the phase difference given in the bottom left. 
The signal in each panel has been normalized to its maximum 
value. 


FIG. 13: Same two-channel sums as Figure 1121 except that 
the events have been summed incoherently. For clarity, the 
prefactor phase information has been marked on panel (e) in 
white. This information can be related to panels 2 and 4 in 
Fig □ The green dotted lines mark nodal lines. The same 
logic can be applied to panel (f), with the difference that the 
4s state does not have angular nodes. 


fringes are comparable to those associated with the “left- 
down” interference. There are no significant diagonal ef¬ 
fects in panel (c) as both distributions are concentrated 
near the axes p n \\ = 0(n = 1,2). This happens because, 
for the two interfering channels, the second electron is ex¬ 
cited to an s state. However, there are some interference 
effects breaking the fourfold momentum symmetry. 

The remaining panels show some implications of chan¬ 
nel interference involving energetically very close and dis¬ 
tant excited states. For the interference of channels 1 and 
5, shown in panel (d), there is a large difference in the 
excited-state energies. This causes small inter-channel 
fringes, so that suppression along the diagonal or anti¬ 
diagonal is not possible. In contrast, for panel (e), the 
excited bound states are energetically very close. This 
implies that fringes stemming from channel interference 
are too thick to cause a diagonal or anti-diagonal sup¬ 
pression. However, the prefactor does this instead (see 
below). For panel (f) the channel-related fringes are even 
thicker, so that the substructure is determined by the 
event interference and the prefactors. 

Figure [13] shows the same coherent sums of channels 


as in Fig. m but, instead, incoherent sums of events. 
All the diagonal and anti-diagonal effects remain and are 
in general stronger without the phase information and 
fringes from different events. The influence of the com¬ 
bined prefactors can also be seen more clearly. For in¬ 
stance, in panels (a) to (c) the features related to p , d 
and s states are very evident, with a further bias intro¬ 
duced by the inter-channel interference. These features 
are (a) probability densities concentrated at the diagonal 
and anti-diagonal; (b) similar probability distributions as 
in (a), but with secondary maxima at the axes; (c) dis¬ 
tributions concentrated mostly at the axes, respectively. 
This happens because the angular momentum quantum 
number l e of the excited states are the same for the two 
interfering channels. 

The situation becomes more complex in the lower pan¬ 
els, in which channels with different angular momenta 
l e are mixed. In this case, we have identified two very 
striking scenarios, which occur for energetically close lev¬ 
els [panels (e) and (f)]. Because the inter-channel fringes 
are very thick in these cases, the shape of the electron- 
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momentum distributions will be mainly determined by 
the prefactors and their phases. For instance, the clear 
diagonal suppression for panel (e) can be explained by 
the phase of the second ionization prefactors, which flips 
at every nodal line. The interfering channels involve ex¬ 
citation to 3 d and 4 p. For a p state, there is one node at 
p 2 1 | = 0, while for a d state there are two. This means 
that, from the top to the bottom of the panel, the phase 
of the p state will flip once, while that of the d state will 
flip twice. Hence, in the second quadrant the channels 
interfere constructively, while in the third quadrant de¬ 
structive interference occurs. The same line of argument 
can be employed for the first and fourth quadrant, but in 
this case the interference pattern will be reversed. This 
shows a case where we have entirely prefactor dependent 
anti-diagonal enhancement. In panel (f), the phase of 
the 4s prefactor will not flip in the momentum region 
of interest, while that of the 3d prefactor will flip twice. 
Hence, this will preserve the fourfold symmetry. Fur¬ 
thermore, interference between the channels will not be 
significant, as d states populate mainly the two diago¬ 
nals and s states lead to distribution localized along the 
axes. This leads to a momentum distribution with peaks 
at the axes and the two diagonals. This distribution is, 
for practical purposes, fourfold symmetric, unless event 
interference is considered [see Fig. H2Tfl]. 

To summarize the inter-channel fringes, prefactors and 
a relative phase can cause a range of interference ef¬ 
fects. The prefactors can cause an inversion, which leads 
to diagonal/anti-diagonal enhancement being swapped. 
They can also apply a phase shift, such that a differ¬ 
ent phase between the two channels is needed for diago¬ 
nal/ anti-diagonal enhancement. Inter-channel interfer¬ 
ence effects are not washed out by more complex superpo¬ 
sitions, as can be seen by directly comparing Figs. fTUIa') 
and (c). 


VI. DISCUSSION AND CONCLUSIONS 

In this work, we have performed an in-depth, semi- 
analytical study of quantum interference in recollision- 
excitation with subsequent ionization (RESI) using the 
strong-field approximation (SFA). Our analysis includes 
interference of symmetry-related features such as elec¬ 
tron indistinguishability, and of different excitation chan¬ 
nels. Overall, we have found that the electron momentum 
distributions are greatly influenced by both types of in¬ 
terference. The main effect of quantum interference is 
to break their fourfold symmetry in Pi\\P 2 \\ plane, while 
the symmetry with regard to the diagonals is retained. 
This fourfold symmetry has been encountered in previous 
RESI studies employing the strong-field approximation 
[25l l26l HU or related methods HU . 

We have shown, by considering a coherent sum of 
symmetry-related events, that interference effects previ¬ 
ously thought to be washed out from integration over 
perpendicular momenta are present in the correlated 


electron-momentum distributions. These effects affect 
the RESI distributions already for a single channel of ex¬ 
citation, via enhancement and suppression near the diag¬ 
onals pi|| = ±P 2 ||- We provide fully analytical expressions 
and estimates for diagonal, anti-diagonal and hyperbolic 
interference patterns. Similar fringes can be seen, but 
have not been explained, in HU , in which RESI has been 
modeled using a strong-field quantum-electrodynamical 
method (see Fig. 7 therein). 

We have also found that inter-channel interference will 
play an important role in RESI, in agreement with the 
results of [32]. We go, however, beyond such studies and 
show that the shape of the electron-momentum distri¬ 
butions will be determined by a complex interplay of 
inter-channel and event interference, and the geometry 
of the excited bound states. This will mainly occur near 
the diagonal and anti-diagonal in the parallel-momentum 
plane. This means that it will mainly affect channels in¬ 
volving excitation to p or d states, while the influence 
on those with s-state excitation will be much less crit¬ 
ical. In particular, for the parameter range of interest, 
the contributions from all channels used in this work are 
comparable. 

We also analyze this interference in more detail using 
two-channel coherent superpositions. In this case, di¬ 
agonal or anti-diagonal enhancement may occur due to 
inter-channel fringes and/or geometry-dependent prefac¬ 
tors. The fringes have hyperbolic shape, and their width 
is inversely proportional to the energy difference between 
the two channels involved. The prefactors will determine 
the region in momentum space to be occupied. In partic¬ 
ular the nodes of the ionization prefactor V P2e of the sec¬ 
ond electron will cause phase shifts, which will influence 
inter-channel interference in specific momentum ranges. 
This enhancement can be manipulated using a relative 
phase. In this context, one should notice that by ap¬ 
propriate choice of channels and relative phase one may 
obtain anti-correlated distributions without resorting to 
bound-state depletion. This latter feature has been em¬ 
ployed in HU in order to suppress the signal in the first 
and third quadrants of the parallel-momentum plane. 

Interestingly, depending on how interference occurs, 
the RESI distributions may exhibit diagonal enhance¬ 
ment (correlation), anti-diagonal enhancement (anti¬ 
correlation), or be spread in the four quadrants of the 
Pi\\p 2 \\ plane. In contrast, for electron-impact ionization, 
the probability density is located only in the first and 
third quadrants and interference effects get washed out 
by transverse momentum integration. This sheds some 
light on experimental findings where different atoms give 
either diagonal or anti-diagonal enhancement (21.]. Di¬ 
agonal enhancement is normally attributed to electron- 
impact ionization. However, for low, below-threshold 
intensities, this could also be related to RESI. Indeed, 
all possible results found experimentally [H, [2(1 HU are 
achievable if we can find the correct superposition of 
channels. Furthermore there exist theoretical studies for 
which anti-correlation has been obtained without excita- 
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tion (38}. This suggests that the ability to manipulate 
diagonal and anti-diagonal enhancement with a phase 
opens up the possibility of control over the RESI pro¬ 
cess, which could lead to various applications. 

Finally, we would like to comment on quantum- 
classical correspondence in RESI. There has been con¬ 
siderable debate in the literature whether NSDI in gen¬ 
eral and RESI in particular is a classical or quantum 
mechanical phenomenon. On the one hand, our results 
show that quantum interference has a striking influence 
on the shapes and localization of the electron-momentum 
distributions. Hence, classic-trajectory methods must 
be viewed with care. On the other hand, highly ex¬ 
cited states may give rise to a quasi-continuum, which 
would allow the existence of a quasi-classical wave packet. 
This would justify the success of classical models. For 
molecules, a larger density of states and enhanced ion¬ 
ization may increase their predictive power [H, [H, [13 ■ 

Furthermore, the SFA considers discrete states and ne¬ 
glects broadening and distortion caused by the field. It 
could well be that these effects lead to a strong overlap 
and thus the creation of a quasi-continuum, washing out 
phase information. However, recent studies of the RESI 
dynamics in phase space have revealed a highly confined 
region that can be associated with trapping in an ex¬ 
cited state [II|]. This would justify using discrete bound 


states and neglecting depletion, and would render inter¬ 
ference important. Another feature which has not been 
included is the influence of the Coulomb potential, which 
changes the topology of the electron trajectories |42l.l43|. 
For a detailed discussion on the advantages and draw¬ 
backs of classical and quantum-mechanical approaches 
in the modeling of RESI see our review article [§j]. The 
present work contributes to this discussion by shedding 
additional light on the role of interference in this context. 
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